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LAKESHORE TWO-DIMENSIONAL DISPERSION 

ABSTRACT 

Hourly two-dimensional dispersion characteristics are deter- 
mined from recording current meter histories for the nearshore areas 
on lakes Erie and Ontario. The current histories were obtained in 
areas within 4 km of shore and at water depths of 10 to 14 metres 
during May to November 1968. A Markov chain process was applied 
to hourly current readings. Three different formulations of the 
stochastic process were tested prior to the selection of the most 
reliable one. The results obtained in applying the developed 
technique compare favourably with results obtained from conventional 
dye injection and drogue studies o 

KEYWORDS: currents, diffusion, instruments and techniques, 
nearshore. 
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LAKE SHORE TWO-DIMENSIONAL DISPERSION 

INTRODUCTION 
Measurements of the nearshore currents in the Great Lakes have 
revealed the highly variable nature of the water movements in both time 
and space. As the currents are the basic mechanism responsible for 
dispersion, it is necessary to understand and develop analytical tech- 
niques to describe the water movementSo The problem can be simplified 
by considering different time and spatial scales separately. This 
technique has been used in air pollution studies (Hilst (4) ), estuarine 
analysis (Okubo (7) ) and Great Lakes Studies (Palmer (9) ). In the 
Great Lakes studies, monthly average dispersion and five hour dispersion 
characteristics were developed from the continuous history of currents 
measured at a fixed point every 10 minutes and averaged to hourly values. 
The five hour diffusion characteristics were determined after lengthy 
analysis on the computer for the four major compass directions^ How- 
ever, the computations are seriously limited to time periods not greater 
than five hours and major compass directions by computer time and 
core space considerations. This represents a serious limitation if 
the technique is to be incorporated as an integral step into a larger 
assimilation model. The purpose of this paper is to outline alternate 
methods for determining hourly two-dimensional dispersion character- 



i sties from current meter records for depths of water 10 to 14 metres 
within 4 km of the shore. These methods require a fraction of the 
computer time compared to the other method (Palmer (9) ). The 
methods also produce maximum, mean and minimum dispersion 
characteristics with the associated probabilities for eight major 
compass directions and various time periods. 

DEVELOPMENT OF METHOD 
From the continuous hourly current records for a month, it is 
possible to obtain a one-hour state probability transition matrix. The 
currents are first classified into 80 designated current states (8 directions 
and 10 speed intervals) as outlined in Table 1. 

TABLE 1 . State number index 



Direction 





337.5° 


22.50 




247.5^ 


292.5" 


Magnitude 


to 
22.5° 


to 
67.5° 




to 
292.5° 


to 
337.5° 


0-2 cm/sec. 


1 


9 




7 


8 


2-4 cm/sec. 


9 


10 




15 


16 


16 - 18 cm/sec 


65 


66 




71 


72 


18 - 20 cm/sec. 


73 


74 




79 


80 



The velocity intervals can be varied to suit the measured velocity 
ranges. An initial state probability vector "q" for the month can 
be obtained. The initial state probability vector "q" is the prob- 
ability of occurrence of each state for the month. A transitional 
state probability matrix can be developed. As an example, suppose 
a loaded die was thrown 1 ,000 times and a record was kept of the 
results. There are six possible outcomes, called states, on each 

throw; namely one, two, three or six. Suppose a two is cast, 

it is possible to determine the probability of all the other possible 
states occurring on the next throw by examining the record to obtain 
what states followed a two state and how often it occurred. A 
different set of probabilities would be obtained by considering a 
different start or initital state; namely one, or three, Considering 
all possibilities produces a transitional state probability matrix six 
by six. The one- hour state probability transition matrix can then be 
obtained from the analysis of the current histories classified into 
states (Palmer (9) ). A sample of the resulting 80 x 80 transition 
matrix appears in Table 2A. 



TABLE 2A. One hour transition probability matrix 



Final State 





SPEED 




to 4 cm/sec. 


2 to 4 cm/sec. 






1, 

to 2 
cnv/sec. 


ANGLE 


337.5° 292.5° 
to to 

22.5° 337.5° 


337.5° 292.5° 

to to 
22.5° 337.5° 


■~ 


0) 
■l-> 
(0 
■t-> 
CO 

•— 1 


337.5-22.5 
292.5-337.5 




A 




(0 

5 


2 to 4 
cm/sec. 


337.5-22.5 
292.5-337.5 




1 








i' 




etc. 





Each block represents an 8 x 8 matrix and a sample representing block 
"A"may be found in Table 2B. 



TABLE 2B. One hour transition probability matrix 







Final State 
2.0 to 4.0 cm/sec. 






337.5° 

to 
22.5° 


22.5° 
to 
67.5° 





292.5° 

to 
337.5° 


Initial State 
to 2 . cm/sec. 


337.5° 

to 
22.5° 


0.01737 


0.00668 





0.00240 


22.5° 

to 
67.5° 


0.00634 


0.01692 





0.00121 


- 


- 


. 





- 


292.5° 
to 

337.5° 


0.00626 


0.00521 




0.01147 



Hence by selecting an initial current a row is specified and the prob- 
ability of any state one hour later may be found by looking at the inter- 
section of that row and the appropriate column. For example, the 
first element of block "A" represents the probability of going from a 
state of 0-2 cn%/sec. and 337.5° to 22.5° to a state of 2-4 cm/sec. 
and 337.5° to 22.5°. 

Three alternate methods were developed using the above 
techniques to predict dispersion characteristics. As a first step, it 
was necessary to develop the state probability vector after a certain 
number of hours given the initial state vector and the continuous 



current meter records. The three approaches Involve different methods 
of determining this final state vector but use the same method of pre- 
dicting the dispersion characteristics. 

The first method involved combining the initial state probability 
vector "q" for the month (essentially a frequency analysis of a month's 
current meter record by states) with various powers of the one hour 
transition probability matrix "T mn^-^^" where i represents the power 
of the matrix (Bharucha-Reid (1); Kemeny (5) ). 

m = 1 

where q = final state probability vector 

after "k" hours, "n" is the 
final state. 

^^ ^ initial state probability vector. 

" m" is the initial state . 



T r«„(l) - one hour state transition prob- 
^ mn ' Yh 

ability matrix raised to "i*-"" 

power where i = k. 
The obvious weakness of this method is that any errors present in the 
first power of the transition matrix will be propagated and intensified 
by the process of raising the matrix to power "i" . 

The second laternative method involves the use of several first 
order transition matrices based on varying the intervals between suc- 
cessive currents used in computing the matrix. For example, a matrix 



could be found hy considering only the probability of transitions 
involving every fifth hourly reading and hence a first order matrix 
for five hours could be found. These matrices could then be com- 
bined with the initial state vector to give the final state vector 
after any given number of hours. 

%(k) = £ Q.T^(k) (2) 

m = 1 

where ^^M - final state probability vector 

after "k" hours, "n" is the 
final state. 

= initial state probability vector. 
" m" is the initial state. 



*m 



T (k) = state probability transition matrix 
after k = 1, 2 , 3 hours. 

This method avoids the necessity of raising the transition matrix to 

higher powers and thus eliminates errors found in the previous method, 

However, computer time increases substantially when it becomes 

necessary to search the monthly records many times to establish the 

various first order matrices. The method is weakened somewhat in 

theory because a first order Markov chain implies that a result is 

dependent only on the directly preceding result. While it is known 

that currents occurring in a particular hour are related to currents 

occurring one hour previously, the relationship of currents separated 

in time by five or eight hours is dependent on the physical size of the 



area as well as the driving force. It is not possible to state as a 
general rule that the current occurring now is dependent on the current 
which occurred five hours earlier. 

The final method tested avoided the weaknesses of the previous 
two methods. The initial state vector and the transition matrix are 
defined as in the first method. However, the higher order final states 
are computed using the final state vector for the previous hour as the 
initial state vector for the present computation. For instance, the 
final state vector for the first hour would be multiplied by the first 
order one- hour transition matrix to give the final state vector for the 
second hour (Gablinger (3) ). 

80 

n •^"» fn rnn 

m = 1 

where q (k) = final state probability 

vector after "k" hours 

Qn^Ck-l) = final state probability 

vector after "k-1" hours 

T (1) = one hour state probability 
transition matrix. 

This method directly relates the results of one time period to subsequent 

time periods. 

The question of how many hours must be considered for short 

term dispersion characteristics can be resolved by determining the 
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number of hours required to approach a monthly mean. In Figures 1 
and 2, the state probabilities are plotted against the various states 
for west and north respectively. The monthly mean state for the west 
direction obtained from a frequency analysis of the monthly record is 
state 23. Figure 1 demonstrates the tendency of the most probable 
state to approach state 23 by the eighth hour, whereas a most probable 
state for the north direction is state 1 which is achieved in the first 
and successive hours. As the various directions at a particular location 
are related, it is necessary to consider the longest time interval 
required to approach mean conditions in any direction. It has been 
found that a period of eight hours is a judicious time interval for 
short term dispersion characteristics. 

A measure of the dispersion chairacteri sties for each hourly 
interval can be obtained by using the state probability for that time 
interval as weighting factors. The weighting factors multiplied by 
the respective states and extended to the time interval considered 
will generate the most probable distances travelled in the eight major 
compass directions for the time period considered. As an example, 
the distance travelled in the south direction after one hour is com- 
puted as follows: 
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State 
(A) 

1 

9 
17 
25 
33 

41 
49 



Mean 






Velocity 






cny'sec. 


Probability 




m 


(C) 


(B) X (C 


i 


0.0132 


0.0132 


3 


0.0081 


0.0162 


5 


0.0057 


0.0285 


7 


0.0034 


0.0238 


f 


0.0009 


0.0081 


11 


0.0002 


0.0021 


13 


0.0001 


0.0013 



Sum = 0.0932 



Most probable distance 
travelled after one hour 



= 0.0932 X 60 X 60 
= 336 cm. 



A comparison of the results by using equations 1, 2 and 3 for 
the mean eight hour pattern appears in Figure 3. It is observed that 
there is little difference in employing the three methods although the 
third method was selected as being the more reliable when considering 
the propagation of errors. A plot of the successive hourly patterns 
produced by the third method appear in Figures 4 and 5 for a location 
on Lake Erie and Lake Ontario respectively. The areas defined by 
the various related hourly contours are representative of the dispersion 
characteristics for that time period. The patterns DO NOT represent 
actual dispersion plumes as the probabilities have been combined 
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with distances travelled in any time period. A two-dimensional 
coefficient " £ " can be obtained by dividing the area contained 
by the related contours by the time period considered. A sample 
of the diffusion coefficients computed by this method appears in 
Table 3 . 



TABLE 3. Diffusion coefficients for Nanticoke 
September 1968 





£ 




Diffusion Coefficients 


cm 


2/sec 




0.67 


X 


104 


1,6 


X 


io4 


3.7 


X 


10^ 


5.7 


X 


104 


8.3 


X 


104 


11.1 


X 


104 


13.9 


X 


104 


17.4 


X 


10^ 



Time 
Hrs. 

1 
2 
3 
4 
S 
6 
7 
8 

It is observed that the diffusion coefficients increase with time. 
This is an expected trend as the longer time periods permit the larger 
scale physical phenomena to become operative and increase the dis- 
persion characteristics, e.g. lake wide wind seiches and littoral 
drift effects become operative at periods of greater than four hours 
on the average. A comparison of the diffusion coefficients obtained 
by this method with results from other studies appear in Table 4. 
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TABLE 4. Results of other studies 



Reference (cmVsec.) 



Csanady (2) 4 x 10^ (one at 2 x 10^) 

Okubo (8) 3 to 6 X lO"* 

Noble (6) 2.44 X 10^ 

The mean conditions for four different locations on lakes Erie 
and Ontario appear in Figure 6. It is observed that a power relation- 
ship exists between the log (time) and the log (area). This is a 
similar situation determined by Okubo (7) in his estuarine study 
where the diffusion coefficient increases with distance from an 
estuary mouth, if one assumes that distance and time are directly 
related. The power varies from .389 to .545 which indicates the 
dependence of diffusion on both the local topography and the time 
of year. 

It is also possible to obtain a measure of the minimum and 
maximum likely dispersion characteristics from the meter records. 
This can be achieved by considering various portions of the state 
probability vector at the end of various time intervals. For instance, 
a measure of the minimum dispersion characteristics could be obtained 
by considering only states 1 to 8. The probabilities for these states 
can be summed and weighting factors can be obtained by increasing 
the probabilities to sum to one. The distances travelled for each time 
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interval can then be obtained by multiplying the weighting factors 
by the state velocity Intervals and the appropriate time intervals 
considered. The probability of this condition is defined by the sum 
of the original probabilities for these states. Similarly, a maximum 
condition can be determined by considering only states above state 41 
and performing a similar operation to determine the weighting factors. 
A selection of the states to be considered for minimum and maximum 
condition will obviously be dependent upon the current history for a 
particular area. A comparison of the results for equations 1 , 2 and 3 
for the minimum and maximum eight hour patterns appears in Figures 
7 and 8, For the minimum there is little difference between the three 
methods, however, equation 2 does generate a significantly different 
pattern for the maximum condition . There are some theoretical problems 
in selecting method two (equation (2) ) and the third method (equation 
(3) ) was selected to determine minimum and maximum as well as the 
mean patterns. The results of applying equations (1), (2) and (3) to 
the data to determine minimum, mean and maximum eight hour dispersion 
patterns appear in Figures 7, 3 and 8. The areas defined by the 
different methods appeeir in Table 5 . 
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FIG. 10 
FRENCHMAN BAY 

JUNE 1968 
EIGHT HOUR DISPERSION 

2 METERS BELOW SURFACE 

MINIMUM, MEAN AND MAXIMUM 

SCALCt r ■ ie,ooo 



70- IS-QL 



FIG. II 
MINIMUM AND MAXIMUM DISPERSION AREAS 
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TABLE 5. Eight hour area comparisons 







Method 






Equation (1) 
X 10'' cm2 


Equation (2) 
X 10'' cm2 


Equation (3) 
X 10^ cm^ 


Minimum 


5,2 


4.5 


4.7 


Mean 


14.0 


10.2 


14.4 


Maximum 


6,600 


5.200 


6,000 



Examples of the maximum and minimum conditions found using 
the third method mentioned above for an eight hour period appear in 
Figures 9 and 10. Once again, the area is defined by the various 
hourly contours which are representative of the dispersion character- 
istics. The areas defined by this method appear in Figure 11 for the 
maximum and minimum power at four different locations on lakes Erie 
and Ontario. It will be noticed that a power relationship exists once 
again for a maximum and minimum condition, but that the slopes do 
not vary as widely as those of the mean, although the intercepts are 
more variable. 



SUMMARY AND CONCLUSIONS 
The Markov chain analytical technique has been demonstrated 
as a reasonable method for determining two-dimensional dispersion 
characteristics from recording current meter data. The figures of 
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patterns DO NOT represent actual plumes. It has been tested in three 
different forms, namely, the development of higher order chains utilizing 
the power of first order chains, higher order chains developed for differ- 
ent time Intervals, and step wise, first order Markov chains applied 
hour after hour. The latter technique has been found to be the most 
reliable to predict hourly two-dimension hourly dispersion patterns 
However, the difference between the three methods is, in most cases, 
negligible. The results obtained by employing the Markov process 
compare favourably with dispersion characteristics determined by 
conventional techniques, such as dye dispersion runs and drogue 
trackings. Further, the Markov chain process can be particularly 
useful in determining the limits of short term (hourly) dispersion 
characteristics. 

A linear relationship exists between the log (time) plotted 
against the log (area) with slopes varying from 0.389 to 0.579 and 
intercepts varying 0.026 to 1 , 1 hours. This is a similar result 
determined for estuarine dispersion. No significant relationship 
could be developed relating diffusion characteristics with month or 
geographic location from this data. The highly variable nature of 
the nearshore area dispersion characteristics of lakes was confirmed. 
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NOTATION 



The following symbols are used in this paper: 

= two-dimensional dispersion coefficient 
cm V sec. 

qj^(k) = final state probability vector after "k" 
hours 



»m 



T mn^*^^ 



initital state probability vector 

state transition probability matrix for "k" 
hour intervals between readings raised to 
the power "i" . 



current 



a horizontal (two-dimensional) movement 
of water which has magnitude (crrv/sec.) 
and direction (degrees) 



diffusion coefficient 



a measure of the mean spread of particles 
transported by water movements or a 
measure of dilution (cmVsec.) 



initial state 
probability vector 



littoral drift 



the probability of occurrence of all possible 
states for a month arranged as a column or 
row vector 

the water movement generated from the inter- 
action of waves and shorelines 



Markov chain 



an array of probabilities which relate the 
state existing now with the state which will 
occur in the next time period. It is a method 
of describing a sequence of related events. 



spatial 



state 



a term relating physical or geometrical 
separation 

physical possible condition that a variable 
can assume 



29 



e.g. 1 C Is a temperature state 

2 cm/sec. at 170° is a current state 
300 to 320 umhos/cm is a conductivity 
state 



state transition 
probability matrix 



similar to Markov chain but restricted to a 
two-dimensional array 



wind seiche 



a periodic water movement resulting from the 
interaction of winds and water surfaces in 
areas with finite dimensions 
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